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Abstract. We present an improved version of Berry's ansatz able to incorporate exactly the existence of 
boundaries and the correct normalization of the eigenfunction into an ensemble of random waves. We then 
reformulate the Random Wave conjecture showing that in its new version it is a statement about the 
universal nature of eigenfunction fluctuations in systems with chaotic classical dynamics. The emergence 
of the universal results requires the use of both semiclassical methods and a new expansion for a very old 
problem in quantum statistical physics. 
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1 Introduction 

During the Chladni meeting, we have honored Chladni's 
curiosity about the complex and beautiful patterns vis- 
ible even in very simple vibrating systems. His interest 
was rewarded not only with the esthetic pleasure we are 
sure he experienced when watching this patterns emerge 
■ by adding some sand grains to his vibrating plates. Quite 
likely he also enjoyed applying a deeper geometrical un- 
derstanding of the modes of vibration (as we would call 
them now) and ultimately the transformation of sounds 
into forms. This curiosity is still with us, partially be- 
cause our modern techniques and tools keep pushing the 
connection between geometry and vibrations into new do- 
mains of science, and partially because some old questions 
remain unsolved. 

During the conference we have seen examples of both 
these connections (ranging from mesoscopic physics to the 
design of musical instruments, or from complex analysis 
to the theory of visual perception) and these old questions 
(such as the ones related with the famous "can you hear 
the shape of a drum?"). 

One connecting link between many of the contributions 
presented at the conference is the subject of this article. 
This connecting concept we refer to is a conjecture and 
a useful practical tool which says that in certain kinds of 
vibrating systems, universality of the corresponding geo- 
metrical patterns appears. Of course, as it stands now, this 
statement is very loose, and a way to formulate it precisely 
is one of the open questions we mentioned above. 



Which systems we are talking about and what does 
universality mean in this context are questions that, thanks 
to the seminal work of M. V. Berry pQ, can be answered 
by now. However, much is still to be done, and we will ask 
(and in a couple of cases answer) other subtle aspects of 
this conjecture. 



2 Berry's conjecture: the random wave model 

In this section we attempt to give a form of the conjecture 
narrowing both the class of systems and the meaning of 
universality. We follow the ideas as presented by Berry in 
his classic paper of 1977 pp. 



2.1 The universal two-point correlation function 

Consider a quantum particle moving inside a 2D domain 
(billiard) fl with area A(f2) and Dirichlet boundary con- 
ditions along its boundary dfi. Berry's main observation 
was that for a given solution ip n (r) of the Schrodinger 
equation 



ipn{r) I r edn = 0, 
the spatial two-point correlation function 
1 
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has a universal form if the classical dynamics inside the 
billiard is chaotic and the points r 1,1*2 are far from the 
boundaries. Here S is a subset of fl with area A(S) con- 
taining many wavelengths. 

The universal two-point correlation function predicted 
by Berry is given by the well known Bessel function pQ 



R n (ri,r 2 ) 



A(f2) 



Jo(k n \ri - r 2 \ 



(3) 



an extremely robust result which has been extensively 
checked both numerically and experimentally [2j. 

There are many arguments (but not a proof) to un- 
derstand this universality. A useful and straightforward 
way is to note that in the semiclassical regime the eigen- 
functions of classically chaotic quantum billiards should 
appear structureless, isotropic and roughly homogeneous 
due to the lack of structure of the classical phase space. 
This observation can be made precise by assuming that 
tpnij") is locally written as a superposition of plane waves 
with fixed wavenumber, directions uniformly distributed 
and random phases. The ensemble 
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with 8 — (Si, ... , Sl) independent random variables uni- 
formily distributed over (0, 27r] and A/" a normalization fac- 
tor is then called a Random Wave Model (RWM). 

It can be easily shown that in the limit L — > 00 the 
RWM correlation function (after average over the random 
phases) is given again by Eq. Q . This is the fundamental 
observation made by Berry pQ , with very deep and power- 
ful implications: the two-point spatial correlation function 
of eigenfunctions associated with classically chaotic dy- 
namics (also called "irregular" ) is universal and given by 
the correlation function of an ensemble of random waves 
with fixed wavenumber. Note, however, that the two kinds 
of ensembles are quite different. On one hand, we have the 
spatial average in Eq. {2]), on the other hand the average 
over random realizations of the phases in Eq. (j4|). In fact, 
for a given billiard shape, there is a third kind of average, 
the spectral average, where the fluctuating variables are 
the values of eigenfunctions with different energy index n. 
This last kind will be of importance later on. 

The importance of the relationship between irregular 
eigenfunctions and random waves goes far beyond this 
analogy concerning the two-point correlation function. The 
point is that the semiclassical analysis used by Berry to 
derive the Bessel correlation result Eq. ^ cannot be ap- 
plied to more general objects. For example, there is no 
known way to relate objects like the fourth-order spatial 
correlation function 



Rn(?i, r 2 ,r 3 ,r 4 ) 



1 



A(S) 



(5) 



ip n (ri + q)il> n (r 2 + q) ip n (ra + q)i> n {ri 4- q)dq 



to the Bessel two-point correlation. The reason lies deep in 
the very construction of the semiclassical approximation 



to the quantum propagator, and it is finally connected 
with the wrong analytical structure of its semiclassical ap- 
proximation. This interesting subject is beyond the scope 
of this contribution and we recommend, e.g, [3] for further 
information. 

After mentioning that the semiclassical techniques can- 
not be directly used to calculate higher order correlations, 
we are ready to formulate Berry's conjecture. Berry's in- 
sight went far beyond the two point case to say that all the 
statistical properties of the eig en junction's spatial fluctua- 
tions in chaotic systems are described by a superposition of 
waves with fixed wavenumber and random phases. In this 
form the Random Wave Model in principle allows for the 
calculation of not only the two-point correlation function 
but any statistical property of irregular eigenfunctions, as 
we show now. 



2.2 Random Wave Model as a Gaussian Random Field 

In principle, in order to calculate higher order correlation 
functions using the RWM, we must insert our definition 
(|4]) into the expression to be averaged (for example Eq. 
(O), then we evaluate the average over the random phases 
and consider the limit L — » 00. This procedure can be 
of use for very simple averages, but it becomes rapidly 
complicated with the order of the statistical measure. 
In general, for a M-point statistical measure 



F^r M (r 1 ;6) 1 ...^r M (rM;S)) 



(6) 



depending on the values of the random wave at fixed 
points i*i, . . . , tmi we must evaluate the integral over the 
random phases 



(F)= I F(i>« WM ( ri ;6),...,t 



F(v)P{v)dv, 



1, RWM 



(r M ;S))dS 



(7) 



where d5 = d8\ . . . d&L ■ We also introduced the vector v 
(v\,..., vm) and its joint probability distribution 



P(v) = lim / 

JV^oo J 



M 



II ^ 



d5, (8) 



in general a very hard calculation. 

Luckily, explicit calculation of the distribution P(v) 
for eigenfunction values is not necessary. Due to a straight- 
forward application of the central limit theorem [2J, the 
distribution of eigenfunction values is Gaussian, namely, 



P(v) = 



1„. T-l 



v /(2 7 r) M det J 



(9) 



where the elements of the so-called correlation matrix J 
(the only set of parameters defining the distribution) are 
given by 

Ji,j = MK\r l - rA). (10) 



Juan Diego Urbina, Klaus Richter: Random Wave Functions with boundary and normalization constraints: 



3 



Equations ([9]) and (fTUj) represent the second (conjectured) 
ingredient making Berry's conjecture such a powerful tool, 
and another source of universal behavior in the statistical 
description of chaotic eigenfunctions. Not only the two- 
point correlation function is universal, also the distribu- 
tion of eigenfunction's values is universally described by a 
Gaussian distribution uniquely fixed by the universal two- 
point correlation function. The third version of the RWM 
conjecture now says that eigenfunctions in quantum sys- 
tems with chaotic classical dynamics have the same statis- 
tical properties as a Gaussian random field with a univer- 
sal correlation function. This is the form of the conjecture 
adopted in more recent works [4]. 

We have seen during this meeting several examples of 
the robustness and power of this conjecture. Not only the 
two-point correlation function has been checked, but also 
the impressive agreement of the RWM with very com- 
plicated functionals of the eigenfunctions such as vortex 
correlations and nodal densities shows the universality of 
the eigenfunction distributions. 

However, as refined calculations near boundaries were 
required (for example in [6]), two major issues moved into 
focus. Both of them demanding changes in the requiere- 
ments to the RWM, and both of them producing devia- 
tions from universality and/or Gaussian statistics. In the 
rest of our contribution we will study these modifications. 



3 Including boundary conditions into the 
Random Wave Model 

Our first concern to describe realistic eigenfunctions be- 
yond the universal isotropic results of the RWM is the 
inclusion of boundaries into the model. Inasmuch as the 
RWM is based on two concepts, namely, a universal two- 
point correlation and an universal probability distribution, 
one has to understand how the boundaries could affect 
each of these ingredients. 

This question was posed already in the frame of the 
statistical description of wavefunctions in disordered sys- 
tems, and the answer appeared to be quite drastic: the 
presence of boundaries is not compatible with Gaussian 
statistics [7]. 

This result has delicate consequences, particularly since 
much of the power and usefulness of the RWM stems from 
the Gaussian statistics which came with a set of power- 
ful rules to perform calculations (most of all Wick's the- 
orem). Is it possible to keep at least by now a Gaussian 
theory but to introduce boundary effects in a less destruc- 
tive way? Contrary to the expectations imported from the 
theory of disordered systems, the answer to this question 
is positive, but we need to introduce more carefully the 
kind of averages our new approach can handle. 

We have seen in section 12.11 that the original RWM 
dealt with spatial averages. It is quite clear that spatial 
averages will in general destroy any information about the 
boundary, unless they are done in a very special way. For 
example, if the spatial region S we use to perform the inte- 
grations in Eq. ^ is circular, we will lose any information 



on the direction where the boundary is located. In fact we 
will reinforce the Bessel result, as shown in [TP] , 

We conclude that any kind of well defined directional 
or non isotropic effect is visible only when the average is 
done in such way that no spatial integration is required. In 
the case of disordered systems this is done by definition, 
since the different realizations of the disordered potential 
provide a natural ensemble, but in the case of a single 
system (as it should be for clean systems) there is no such 
mechanism. 

The only option left is to no longer deal with a single 
eigenfunction and perform the averages for fixed sets of 
positions, without any spatial integration and varying the 
energy index n. As we will see, this idea will open the door 
for a very compact and appealing formulation of Berry's 
conjecture. 



3.1 Eigenfunction's averages over spectral windows 

Consider now a functional F]%1>] depending on the set of 
M positions r±, . . . , tm as parameters through the values 
of the function at those points 

Fty] =F(^(r 1 ),...,^(r M )). (11) 

When we plug different eigenfunctions into F, we obtain a 
fluctuating quantity. In a natural way we define the spec- 
tral average of the functional F around energy e, indicated 
by calligraphic symbols from now on, as [5] 

•F=— ^5>(e-e n )F[^ n ] (12) 

where w(x) is a normalized window function around x = 
and p w (e) is the density of states smoothed over a window 
w. Our objective is to construct a random wave function 
to describe such averages, inasmuch as Berry's RWM de- 
scribes spatial averages, and based on similar statistical 
assumptions. 

First we note that our spectral average indeed respects 
boundary conditions. For example, if one of the positions 
in the set r\, . . . ,Tm lies exactly at the boundary, then 
all the eigenfunctions must vanish there exactly, not only 
on average but sharply. How can we implement such hard 
constraint into the otherwise very smooth Gaussian distri- 
bution? The answer is actually quite simple: do not change 
anything. 



3.2 An ensemble of random waves satisfying arbitrary 
boundary conditions 

Consider the exact two-point correlation function associ- 
ated with the spectral average in Eq. (|12p. 

R w (r 1 ,r 2 ,e) = — -r^ ^2 w ( e ~ e n )^n(r 1 )^ n (r 2 ), (13) 
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and note that it satisfies exactly the boundary conditions. 
Now, let us construct the following Gaussian distribution 



Pr(v) 



1 



-±l>.RT 



V(27r) M det R 



fixed by the system-dependent correlation matrix 

RiJ Rwy^ii^ji^ji 



(14) 



(15) 



depending on the set of positions T\, . . . , Tm and the en- 
ergy e as parameters. 

It is not difficult to check that the probability distri- 
bution Pr(v) converges to a delta function 8{vi) when the 
point Ti lies at the boundary. Namely, Eqs. ([HI [TS|) de- 
fine an ensemble of random waves satisfying sharply the 
boundary conditions. Wc mention that, although we con- 
sidered the Dirichlet case for simplicity, any kind of Hcr- 
mitian boundary condition can be treated exactly in the 
same way. 

Using our construction, we can present yet another 
version of Berry's conjecture, now with the extra gain in- 
cluding arbitrary boundary conditions: spectral averages 
of junctionals in classically chaotic quantum systems are 
universally described by a Gaussian distribution fixed by 
the system- dependent two point correlation function. 

At this point, a natural question arises, namely, what 
is the connection between the universal predictions of the 
RWM and the system-dependent predictions of our new 
ensemble? Luckily, they are compatible, since standard 
techniques based on the semiclassical expansion of the 
propagator can be used to calculate the correlation func- 
tion, and R w (ri,rj,e) can be represented in a multiple 
reflection expansion of the form [§UTD] 



Rw 7 Tj ; ^) 



1 



-Jo{k(e)\n - Tj\) 



sum over reflections at the boundary 



(16) 



Ri.b(ri,r j: e) 
LiLj sec 9 



1 



J Q (k(e)\ri 



2k- 



A(Q) 

1/2 cos [k(e) (Li + Lj)- 



(17) 



Li + La 



A(f2)^2ir(L l + L 3 ) 



representing the effect of a single reflection with the bound- 
ary at point r with incidence/reflection angle (9, we can 
explicitly calculate the effect of such boundary on any 




Fig. 1. Exact calculation (symbols), universal result (horizon- 
tal line) and one-bounce approximation (continuous line) of the 
average intensity {ip 2 (f2)) = ^?iu(j"2, f 2, e) in the cubic confor- 
mal billiard [12]. The average intensity is plotted as a function 
of the distance r | = \r? — r i to the origin ri while moving the 
observation point r2 along the line indicated in the inset. Non 
universal effects are clearly visible and well described by the 
one-bounce correlation function of Eq. (|17[1 (diffraction effects 
must be incorporated). Note in particular how the boundary 
condition (Dirichlet) is respected 



Therefore a similar expansion of arbitrary averages around 
the universal RWM results can be constructed. This method 
was used in [8] to formally show the equivalence between 
the theory of eigenfunction statistics in disordered metals 
and the improved form of Berry's conjecture. In particu- 
lar, all the results of the constraint Random Wave Mod- 
els proposed to incorporate simple geometries (as in [llj l 
are recovered, extended and generalized using the random 
wave model with boundaries [i"3] . 

For example, using the Gaussian distribution Eq. (| 14[) 
together with the following ( "one bounce" ) two-point cor- 
relation function 



statistical measure 



k and L 



iti) 



10] . Here the local curvature at r is 
— 7*1. The extra phase <f> is related to 



the topology of the local classical flux around the reflected 
path and can be also explicitly constructed [14]. 



Figure |T]) shows the average intensity (ip 2 (r2)) (given by 
Rw(i f '2, i"2, e) in Eq. (|13pl as a function of the distance 
|r| = \r2 — t*i| from the origin r*i to the observation 
point r2 as it moves towards the boundary of a partic- 
ular chaotic billiard, shown in the inset. Symbols are ex- 
act numerical results obtained from the numerically ex- 
act eigenfunctions, while the continuous line is given by 
Ri.b(f2,r2,e) (including diffraction effects). As it can be 
seen, the presence of the boundary drastically modifies 
the local structure of the averages when compared with 
the 1/A(f2) ~ 1/tt constant result. 

Our ensemble of random waves with boundaries used 
our freedom to incorporate non-universal effects by means 
of the two-point correlation function, leaving as the only 
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conjecture a statement about the universal Gaussian fluc- 
tuations of irregular eigenfunctions. However, to construct 
an ensemble respecting all the properties of the exact 
eigenfunctions, every member of the ensemble should be 
properly normalized. As we will see, Gaussian statistics 
is not compatible with this constraint, and the Gaussian 
conjecture must be modified. 



4 Projecting the ensemble of random waves 
onto the set of normalized functions 

Our goal in this section is to construct an ensemble of 
random functions, respecting sharply the boundary con- 
ditions, such that every member of the ensemble is ex- 
actly normalized. This is more than an academic question: 
it is fundamental to have normalization correctly taken 
into account, because otherwise the results of the ensem- 
ble average will contain spurious fluctuations affecting any 
statistic beyond M = 2, as we shall explain. 

The so-called normalization problem was first noted in 
[5] where the corrections due to normalization for disor- 
dered metals were found. In the last section of this con- 
tribution, we will show how to obtain such corrections in 
the appropriate limit of our re-formulated Random Wave 
Model, giving strong support to our construction. 

In order to appreciate the kind of spurious effects which 
we get using Gaussian statistics, we consider the func- 
tional 



F[r/>] = J *{j 2 (r)dr, 



(18) 



which should have zero variance if every member of the 
ensemble is properly normalized, since we have F[ip] = 1 
over the set of normalized eigenfunctions. However, let 
us calculate explicitly this variance by using our random 
wave model: 



(F 2 ) - (F) 2 = J (ip 2 (r)ip 2 (r'))drdr' - 1 (19) 
= 2 J R 2 w (r,r',e)drdr', 

where we have used the identity J R w (r, r, e)dr = 1. Re- 
sult Eq. (fT9]) is obviously not zero, and this indicates that 
if we calculate the spectral average of a functional of order 
higher than two, fluctuations of the normalization integral 
will produce spurious contributions. 

The mathematical problem is the search for a prob- 
ability distribution over the set of normalized wavefunc- 
tions, as close as possible to the Gaussian (since we have 
very strong evidence for it), and uniquely fixed by the 
given two-point correlation function (since this object is 
calculated exactly and does not suffer from normalization 
problems). Of course "as close as possible" should be de- 
fined properly, but presently the structure and topology 
of the space of distributions living on the set of normal- 
ized functions and obeying a fixed two-point correlator is 
still poorly understood; so following an educated guess, 
we will try to construct one of such distributions hoping 



that its predictions match the results for disordered sys- 
tems. However, before carrying out this program, we need 
to take a further step. 



4.1 The non-local and non-stationary character of the 
normalization constraint 

One of the convenient properties of the Gaussian distribu- 
tion Eq. (|14|) is its local structure. Locality in this context 
means the following: let us define 



v N 

VM 



,ip(rN), ■ 



with N < M, and consider a functional G[ip] depending 
on the subset of positions n, . . . , r^. It is then true that 
its Gaussian average Q satisfies 



Q = 



1 



l=/ GW e-W^ Mrf „ M(2 0) 
G(v N )e-i*> 



v /(2 7 r) Ar det R' 



R' .v„ dv 



where R' is a correlation matrix obtained by simply re- 
moving the last M — N indeces from R (from now on we 
suppress the w index of the correlation matrix). In other 
words, the average of a functional depending on a given 
subset of positions in still Gaussian (stationarity) and de- 
pends only on the eigenfunction correlations among only 
those points (locality). 

The locality and stationarity of the Gaussian distribu- 
tions are at the very heart of its generality. They mean 
that the distribution has the same form independent of 
how many points the functional has and that the averages 
depend only on the points where the functional is defined. 
Because of this, the most general form of the Gaussian 
probability distribution is the continuous limit of Eq. (TT4|) , 
that is, a functional distribution 



PrUPXDUP] = lim 1 = 

11 1 1 M->™ ^(27r) M dct R 



(21) 



x e "i YZ ^^"^^^(n) . . . #(r M ) 
Mm 



(27r) M det R 



where, after proper scaling of 4>( r i) with M absorbed by 
A/m, the quadratic form in the exponent is 

M 

(V-lirV) = Jim -^(nXR- 1 )^,) 

M— >oo 1V1 * — ' 

id 

V<r)(R -1 )(r,r')V>(r')dr<i!r', (22) 



while R (referred from now on as the "covariance" or 
"correlation" operator) with matrix elements (R)j* = 
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(ri\R\rj) is naturally denned through the spectral win- 
dow function w(x) and the system Hamiltonian H as 



R 



w(e - H) 
Tr w(e — H) 



(23) 



The notations [0]ij, (P)(n,rj) and (ri\0\rj) will be used 
from now on to indicate matrix elements. 

Our expressions (|21j22|) and (|23|) have the very impor- 
tant property of being independent of the representation, 
in the sense that they define the probability of finding a 
given function ip, not a particular set of values for their 
projection along a particular basis. 

Now let us consider the effect of normalization con- 
straints on the probability distribution in Eq. (|2~T1) . The 
requirement that all the members of our ensemble are nor- 
malized necessarily implies that the correct distribution 
that we indicate with Pjj AP (anticipating its "Gaussian- 
Projected" character) must be sharply concentrated on 
the unit sphere in the space of functions, 



pGAP 



WPW\ ~ - i)v[ 



(24) 



From this condition we see that, whatever we do, the nor- 
malization constraint destroys the stationary and local na- 
ture of the distribution, since any set of points is connected 
with any other through the normalization. In other words, 
and in sharp contrast to the Gaussian distribution, aver- 
ages of functionals defined over any subset of positions will 
feel the fluctuations of the eigenfunctions at any point of 
the space. The ensemble of normalized random functions 
cannot be neither local, nor stationary. 

However, evidence shows that a Gaussian (local and 
stationary) distribution must be a very good approxima- 
tion, and we will see that we can indeed construct a dis- 
tribution over the unit sphere as close to the Gaussian as 
we wish. Interestingly, although we start talking about a 
single particle inside a quantum billiard, the solution will 
come from the field of quantum statistical physics. We 
mention that another intersection between random waves 
and quantum statistical physics is presented in [15] . 



4.2 The connection with quantum statistical physics 

Besides being defined over the unit sphere in the space of 
functions, our distribution must fulfill the basic require- 
ment that the associated two-point correlation function is 
known: 



^(riMr^Pg^MPM = (n\R\ r] ). (25) 



Equations (|24I25|) define our problem. Not much is known 
about the possible solution, but at least it is easy to see 
that it is not unique [16j . At this point our hope is that 
the physics of irregular eigenfunctions can help with extra 
conditions when looking for a unique solution, but this 
issue remains open. 



Surprisingly, exactly the same question has had a long 
history in the field of quantum statistical physics (we fol- 
low [16 ). In this context, one is interested on the possi- 
ble distribution P p [tf>] of individual normalized wavefunc- 
tions for systems in thermal equilibrium characterized by 
a given density operator 



-0H 



Tr e-Pi* 



(26) 



where H is the Hamiltonian and the inverse temper- 
ature. Consistency with the results of quantum statisti- 
cal physics demands that averages of expectation values 
(^|0|^i), weighted by P p [i(>], should give the usual result 



(iP\0\^)P p [i;]V[^] = Tr (pO), 



for arbitrary O. This in turn implies 



\P\rj), 



(27) 



(28) 



and then it is clear that P p and P^ AP satisfy the same 
conditions, but in the first case we talk of a "correlation 
matrix" instead of a "density operator". In fact, R is also 
a positive, analytical function of the Hamiltonian with 
unit trace, so it is a density operator corresponding to the 
microcanonical distribution. 

One could think, inasmuch as in the case of few de- 
grees of freedom, that we cannot choose among all the 
possible solutions of conditions (|24I25[) . However, this is 
not the case. There is a guiding principle frequently used 
and well known in the framework of classical and quantum 
statistical physics, known as "Principle of Maximum En- 
tropy" (PME) saying that, for a given set of constraints, 
the best option for a probability distribution is the one 
which maximizes the associated entropy functional. The 
exact formulation and the implications of this principle go 
beyond the scope of this paper, and the reader is encour- 
aged to read the excellent literature available [17] • How- 
ever, as it is shown in [16] . when applied to the present 
problem the PME selects one particular kind of distribu- 
tion, known as the "Gaussian Projected Ensemble" and 
given by 



P 



GAP 



m = / s u 



<0|$Pfl[<flZ>[0. (29) 



We mention that in the microcanonical case (our main 
interest here since it defines the Random wave Model), 
this is the unique solution to the maximization problem. 

Following a long and successful tradition, we assume 
that the methods and techniques of statistical physics can 
be applied, or at least can shade some light, to the na- 
ture of systems with complex behavior but few degrees of 
freedom. We then propose the following version of Berry's 
conjecture, incorporating exactly both, boundary and nor- 
malization constraints: In systems with classically chaotic 



Juan Diego Urbina, Klaus Richter: Random Wave Functions with boundary and normalization constraints: 



7 



dynamics, spectral averages of functionals defined over the 
set of eigenf unctions, are given by the corresponding av- 
erage over the Gaussian Projected Ensemble with fixed 
system-dependent covariance matrix. In other words, we 
conjecture that the statistical fluctuations of irregular eigen- 
functions are universally described by the Gaussian Pro- 
jected Ensemble. 

Of course, only the successful application of the conjec- 
ture will decide whether it expresses or not a fundamental 
property of irregular eigenfunctions, and this is work is in 
progress. There are, however, some predictions that can be 
checked and contrasted with existing results coming from 
the theory of disordered systems (a connection which has 
been very useful in the past) [2]. Before that, however, 
we must be able to operate with the GAP distribution, 
something that has not been done until now due to its 
non-local and non-stationary structure. 



4.3 Correlation functions for the GAP ensemble 

In this section we briefly discuss recent work on the cor- 
relation functions associated with the GAP measure [18] . 
Some of this results are currently used for specific calcu- 
lations, hoping that they can be tested against numerical 
or experimental evidence, and others can be directly com- 
pared with known results for statistics of eigenfunctions in 
disordered systems. To be more precise, we first introduce 
some basic concepts. 

A local functional over the N- dimensional subspace 
Hn of the Hilbert space Tt is a function of the values of 
the wavefunction at the given (finite) set of N positions 



F[ l p] = F^(r 1 ),...,^(r N )). 



(30) 



Using the definition Eq. (1291) . the GAP average of F 
is given by the functional integral 



Fgap — 



VdetA 
M 

VdetR 



e 2 



(31) 



where <, > is the inner product on H, \ip\ — Ktpltp)] 1 ^ 2 the 
norm and J\f is an (infinite) normalization constant. The 
operator (also known as /5 _1 ) is well defined because 
of the properties of the function w(z), which is assumed 
to be analytic and positive 

For future reference, the result of the Gaussian (non- 
projected) average is given by 



VdetR 
F(Vi 



e -iv-(fl«r 1 .v 



which is calculated by means of the usual iV-dimensional 
Gaussian integral over the vector variable ip = (^i, . . . , ipisr) 
The correlation matrix is the projection of the opera- 
tor R on Hat with entries given by 



\R 



N\i,j 



{ri\R\rj) for r^Tj e H N , 



(33) 



while det n is the determinant otiTLm- 

It is essential that for an arbitrary operator A on TL 
and an analytical function f(z), 



f(A) N + f(A N ) 



(34) 



although both sides of the equation are defined over the 
same space Hn- In particular, the matrix (i?/v) -1 is con- 
structed by projecting R and then inverting, rather than 
projecting R . 

The pure technical problem with the GAP averages is 
that, contrary to the Gaussian case, one cannot directly 
integrate out the components of tp belonging to the sub- 
space H./H.N. This problem finally boils down to the ap- 
pearance of the norm \tp\, with components over the whole 
TL, as a denominator in J- gap, as it is seen in Eq. (|3ip . 
However, if F is a homogeneous functional 

F< 2 ">(A^(ri), . . . , \iP(r N )) = A^F^^n), . . JM) 

(35) 

these denominators can be collected and written as an 
exponential by means of the following identity 

1 1 f°° 

W = W)L < 36 > 

valid for Re(^) > 0, where F(x) is the gamma function. 
The integrals in Eq. (|3ip are now purely Gaussian, and 
carefully considering Eq. (j3"4"|) the desired reduction of the 
integration over 7i to integration over Hn is achieved by 
straightforward discretization of the functional integral. 

The results depend on the order of homogeneity 2fi. 
For fi = and fi — 1 we have 



■^GAP ~ ■^G ) (^ N ) ~ '^'da'^G^ + a {R 2 )N)\a=0i 



7(°)/ 



d 



7(0), 



GAP — J G V-"-JV (37) 

The result for 2/i = 2 implies as particular case that both 
the Gaussian and GAP averages give the same result for 
the associated correlation (density) operator, but it is far 
more general. For < fi < 1 and 1 < // < 2 we have not 
been able to calculate the corresponding integrals. 

For fi > 2 we first define a one-parameter deformation 
of the correlation operator 



Q(x) = R(l + 2XR)- 1 , 



(38) 



'N 



(2w) N det N R N 



dxp (32) U2/J,) _ 



where the non-negativity of R assures existence of the in- 
verse for x > and renders Q(x)/TtQ(x) a density opera- 
tor as well. With this definition the GAP average is given 
by the following one-parameter integral 



1 



GAP 



F( M -1) 



J-g(Rn), 



det(l + 2xR) 



^\Q(x)N)dx. 



(39) 
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Equation (|39|) together with the definition Eq. ([38]) is the 
main result of this analysis. It allows us to calculate the 
GAP averages in terms of the Gaussian ones. 

It is important to remark again that normalization ef- 
fects are non-local in the sense that any average respect- 
ing normalization (as the GAP), cannot be defined only 
in terms of the projected density matrix Rn acting on 
Hn- The degrees of freedom associated with the comple- 
ment H/Hn appear through the terms (R 2 )n, {R 3 )n, ■ ■ ■ 
in the Taylor expansion of Q(x). As an example consider, 
for \n), \rj) S H N , 



(r t \R 2 \r 3 ) 
1 



(40) 

ri|^|r>(r|B|rj)dr 



A{(2) 

J R{r u r)R{r,r 3 )dr 



A{(2) 

Eq. (|40[) shows that any dependence non-linear in R will 
connect the points where the functional is defined with 
points in the whole space. Since any GAP average (ex- 
cept the 2/i = 2 case) is a non-linear function of R, and 
not of (R) jv , non-local effects are a generic feature of this 
distribution. 



5 The relation between the Gaussian 
Projected and the Gaussian ensembles 

Now we proceed with the asymptotic expansion of (|39[ 
The key point is the analytical structure of the function 



p(z) = det(l + 2zR). 



(41) 



We note that p(z) is the spectral determinant of the op- 
erator R and then an entire function with zeros at, and 
only at, 

1 



2R( ej 



ej G spectrum of H. 



(42) 



Since the operator R is a non-negative function of the 
Hamiltonian, zj is real negative for all j. This indicates 
that, contrary to x < 0, for x > the function p(x) cannot 
be oscillatory, and the following expansion is uniform in x 



1 



det(l + 2xR) 



exp 



1 



Tr log(l + 2xR) 



(43) 



exp 



{-2x) h 



2 ^ 

k=2 



Tr R K 



2xR) 



Eq. (|4*3"j) is an expansion of the function 1/y det(l 
in powers of x, well behaved upon integration. 

The operation that makes our results divergent is the 
following expansion of the deformed density matrix 



R 



1 + 2xR , 



J2(-2x) k (R k 



k+l\ 



IN- 



(44) 



since its radius of convergence demands 2x < R(ei)^ 1 
while the integration demands x — * oo, rendering the ex- 
pansion of the full integral just asymptotic. 

Equations (|4*3"|) and (|4~4"|) provide a systematic expan- 
sion of the integrand in Eq. (|3"5|) in powers of ir. After term- 
by-term (trivial) integration, we get the desired asymp- 
totic expansion of the GAP average around the Gaussian 
result. It is then of importance to identify the small pa- 
rameters of this expansion, in order to have control over 
the successive approximations for specific situations. 



5.1 Physical interpretation of the small parameters 

In order to identify the small parameter of the expansion, 
we remark again that R is adensity operator and then it 
is non- negative and of unit trace Tr R — 1. Then 



Tr R n < Tr R m , for n > m, 



(45) 



holds, and the asymptotic expansion can be organized 
in ascending powers of the covariance (density) operator, 
namely, the order of a term is given by the total power of 
R. These considerations do not hold for the projected den- 
sity matrix: Rn can appear to any power without affecting 
the order of the corresponding term in the expansion. For 
example, given |r,-), \rj) eHn, 



R 



R 2 
R 



1,3 



is of second order, but 



is still of first order. 



N 



k=0 



In order to interpret physically the small parameter, 
we consider two examples. 

— For the microcanonical case of main interest for us, the 
function w(z) in (|23j) is a window function of width e. 
It is easy to show that Tr R 2 is proportional to A/e 
where A is the mean level spacing. Then we see that 
for the microcanonical ensemble the small parameter of 
the asymptotic expansion is l/L where L is the number 
of energy levels inside the energy window. 

— For the canonical ensemble at inverse temperature fi = 
1/ksT the covariance operator reads w{z) — %~" z , and 
a simple calculation gives Tri? 2 ~ A/ksT, so for the 
canonical ensemble, the small parameter is the ratio 
between the mean level spacing and the thermal energy. 

Having now control over the asymptotic expansion, we 
check whether its zero-order term actually gives the Gaus- 
sian results. 



5.2 Recovering the Gaussian limit 

From the theoretical point of view, the most important 
property of the asymptotic expansion is that it provides a 
way to calculate perturbatively the GAP averages starting 
with the Gaussian case and keeping the approximation 
under control. To see how we obtain to first order the 
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Gaussian average, we truncate both expansions (|43|) and 
(1441) to first order in R to obtain 



1 



det(l + 2xR) 
R 



1 + 2xR 



c + 0(Tr R 2 ) 
R N + 0{R 2 ) 



(46) 



N 



Substitution in Eq. (f3"9")l gives 

x"- 2 e- x ^\R N )dx + 0((R 2 ) N ) 



T 



(2/*) 



1 



GAP 



r( M -i) 

T^\R N ) + 0{{R 2 ) N ). 



Where the definition of the gamma function was used. 
Our result Eq. (gT]) shows that for Tr R 2 -> (which 
means high number of levels inside the window function 
in the microcanonical case or high temperature limit in 
the canonical), the GAP and Gaussian averages coincide. 
Although it looks intuitively trivial, it is of importance 
since it is obtained keeping the approximation under strict 
control. 

All results based on the asymptotic expansion have a 
limited accuracy. Poincare criteria for asymptotic series 
demand those to be truncated just when the typical term 
stops decreasing with the order. It is easy to see that, the 
higher the order 2\i of the functional, the earlier one must 
truncate the expansion. The method fails if not even the 
very first term can be used. This happens for functionals of 
certain order 2/i inversely related to Tr R 2 . This actually 
means that the most interesting effects, those beyond a 
perturbation expansion around the Gaussian result, occur 
for higher correlation functions. 

However, already for small-order correlations truncated 
at the first term of the asymptotic expansion, we obtain 
interesting results, as we will see now. 



R 




R 




l + 2xR 




1 


f 2xR 





(49) 



Next, following the idea of an asymptotic expansion in 
terms of 1/L with L the number of energy levels inside 
the energy window, we use the expansions (jl3")l and (|4"4")) 
up to R 2 ~ 1/L. We get, after substitution in Eq. ([39l 
and calculating the elementary integrals, 



y GA p = y G {R2) 

- 2([R 2 ] iyi [R]jj + [R]j,j[R]i 
+ 0(R 3 ), 



(50) 

[R]^[R} 3J Ty R 2 ) 



(47) where ^(^2) is the usual Gaussian result 

y G {R 2 ) = [R}i4R] j j+2{R} i!j [R] iij . 



6 An example: asymptotic expansion of the 
intensity correlator 

As a non-trivial example of the method, we calculate the 
first correction due to normalization to the intensity-intensity 
correlation function, defined by the local homogeneous 
functional 

Y\i>] = ^{r^irj). (48) 

In this case N = 2, the subspace of interest Hn = is 
expanded by the states {|rj), \rj)} and 2/i = 4. 

In order to apply our ideas, we first identify the Gaus- 
sian result y G (Q(x)N) in terms of the deformed density 
operator Q{x)n- This is a standard result founded by us- 
ing the Wick contractions as 

y G {Q{x) 2 ) = [Q(x)] iti [Q(x)] jtj +2[Q(x)] i<j [Q(x)] iij 



(51) 



It is easy to see that the extra terms in Eq. (1501) render 
the variance of the normalization integral Eq. (|19|) exactly 
zero. This fact appear more clearly when we consider the 
universal limit of the theory, as we do next. 



6.1 Universal limit of the intensity correlator: 
connection with the theory of disordered systems 

Our results so far are the asymptotic expansions of the ex- 
act GAP averages up to second order in the small parame- 
ter. We must remember, however, that for the microcanon- 
ical case our results can be further expanded around the 
universal covariance matrix with entries defined through 
the Bessel function. The interplay of these two expansions 
is a delicate matter, but at least we can see what happens 
when we restrict ourselves to the isotropic result for the 
correlation function. 

In order to find the universal limit of our result Eq. 
(f50|) . we use the Bessel form of the correlation function 
Eq. ([3]), the expression Eq. (|40|) for the matrix elements 
[R 2 ]i.i = (R 2 )(ri,ri) and the continuum version of the 
trace to obtain: 



A{f2) 2 y GAP ~ l + 2J 2 (fc(e)| 



(52) 



A{Q) 
2 

A(Q) 
2 



J 2 {k(e)\ ri ~r\)dr 
J 2 (k(e)\r' -r 3 \)dr> 



J 2 (fc(e)|r -r'\)drdr', 



R 






R 


l + 2xR 




1 - 


V2xR 



J 3,3 



with fc(e) the wavenumber corresponding to the energy at 
the center of the spectral window. 

Eq. (|52p , obtained using the Principle of Maximum En- 
tropy and semiclassical considerations, is in fact the same 
result obtained using the supersymmetry technique in dis- 
ordered systems [5]. The first line represents the usual 
Gaussian result, and the non-local terms, involving inte- 
grations over the whole space, are the normalization cor- 
rections up to second order in the inverse number of energy 
levels inside the spectral window. 
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We note that the theory of wavefunction's statistics in 
disordered systems, known as the diffusive sigma model, 
has a very different mathematical structure and physical 
motivation than our original idea of describing irregular 
eigenfunctions in single quantum systems. However, find- 
ing a way to use Berry's conjecture and/or the Princi- 
ple of Maximum Entropy to derive a result coming from 
the theory of disordered systems points towards a deep 
connection beyond the apparent differences between those 
fields. Research towards understanding the guiding princi- 
ple behind this connection has been done for the universal 
fluctuations in the spectrum of statistical, chaotic and dis- 
ordered systems [2], but not at the level of eigenfunctions. 
The present contribution goes into this direction, guided 
by Chladni's broad and elegant way of thinking. 



17. See, for example, W. Grandy Foundations of Statistical 
Mechanics (D. Reidel Publishing Company, Holland, 1993). 

18. J. D. Urbina, in preparation. 
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